Computational combinatorial ligand design: Application to human α-thrombin

Summary A new method is presented for computer-aided ligand design by combinatorial selection of fragments that bind favorably to a macromolecular target of known three-dimensional structure. Firstly, the multiple-copy simultaneous-search procedure (MCSS) is used to exhaustively search for optimal positions and orientations of functional groups on the surface of the macromolecule (enzyme or receptor fragment). The MCSS minima are then sorted according to an approximated binding free energy, whose solvation component is expressed as a sum of separate electrostatic and nonpolar contributions. The electrostatic solvation energy is calculated by the numerical solution of the linearized Poisson-Boltzmann equation, while the nonpolar contribution to the binding free energy is assumed to be proportional to the loss in solvent-accessible surface area. The program developed for computational combinatorial ligand design (CCLD) allows the fast and automatic generation of a multitude of highly diverse compounds, by connecting in a combinatorial fashion the functional groups in their minimized positions. The fragments are linked as two atoms may be either fused, or connected by a covalent bond or a small linker unit. To avoid the combinatorial explosion problem, pruning of the growing ligand is performed according to the average value of the approximated binding free energy of its fragments. The method is illustrated here by constructing candidate ligands for the active site of human -thrombin. The MCSS minima with favorable binding free energy reproduce the interaction patterns of known inhibitors. Starting from these fragments, CCLD generates a set of compounds that are closely related to high-affinity thrombin inhibitors. In addition, putative ligands with novel binding motifs are suggested. Probable implications of the MCSS-CCLD approach for the evolving scenario of drug discovery are discussed.     

Factors Contributing to Solubility Synergism of Some Basic Drugs with β-Cyclodextrin in Ternary Molecular Complexes

Ternary complexes exploiting solubility synergism (SSn) between basic drugs and β-cyclodextrin (β-CD) in the presence of an organic hydoxy acid have been reported to provide the pharmaceutical technology with highly soluble ternary complexes, even with the least soluble β-CD. In this work, phase solubility techniques were used to study factors affecting SSn in aqueous solution, which may help in understanding the mechanism involved in ternary complex formation in solution, under equilibrium conditions. The equilibrium solubility of both β-CD and each of 8 structurally unrelated drugs were measured in tandem in the presence of different acid types at low and high pHs, and at different time intervals over a period of 1–40 days. The results indicate that SSn is evident regardless of acid type (organic and inorganic) at low pH, but the extent of SSn is acid type dependant and is limited by the drug salt solubility product constant (pK _sp). Among different drugs, no apparent trend exists between drug salt solubility and the extent of SSn, but lowering drug salt solubility by increasing pH depresses SSn. The results also reveal no apparent trend between the magnitude of the complex formation constant (K _ij) and SSn. For example, drugs of low K _ij values such as astemizole, cisapride and sildenafil do not show any SSn, yet ketotifen and pizotifen, which also have low K _ij values, exhibit substantial SSn. However, the solublizing power of β-CD represented by the slope of phase solubility diagram can be used as a marker for SSn (slopes exceeding 0.4 induce SSn).     

Functionality map analysis of the active site cleft of human thrombin

Summary The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors. Abbreviations: MCSS, multiple copy simultaneous search; PPACK, d-phenylalanyl-l-propyl-l-arginine chloromethane; NAPAP, N -(2-naphthylsulfonylglycyl)-d-para-amidinophenylalanylpiperidine; argatroban, (2R,4R)-4-methyl-1-[N -(3-methyl-1,2,3,4-tetrahydro-8-quinolinylsulfonyl)-l-arginyl]-2-piperidine carboxylic acid; rms, root mean square. The thrombin residues are numbered according to the chymotrypsin-based numbering by Bode et al. [8]. P1, P2, P3, etc., denote the peptide inhibitor residues on the amino-terminal side of the scissile peptide bond, and S1, S2, S3, etc., the corresponding subsites of thrombin     

Thermodynamic Parameters for the Association of Fluorinated Benzenesulfonamides with Bovine Carbonic Anhydrase II

This paper describes a calorimetric study of the association of a series of seven fluorinated benzenesulfonamide ligands (C₆H_nF_5?nSO₂NH₂) with bovine carbonic anhydrase II (BCA). Quantitative structure–activity relationships between the free energy, enthalpy, and entropy of binding and pK_a and log P of the ligands allowed the evaluation of the thermodynamic parameters in terms of the two independent effects of fluorination on the ligand: its electrostatic potential and its hydrophobicity. The parameters were partitioned to the three different structural interactions between the ligand and BCA: the Zn^II cofactor–sulfonamide bond (≈65 % of the free energy of binding), the hydrogen bonds between the ligand and BCA (≈10 %), and the contacts between the phenyl ring of the ligand and BCA (≈25 %). Calorimetry revealed that all of the ligands studied bind in a 1:1 stoichiometry with BCA; this result was confirmed by ¹⁹F NMR spectroscopy and X‐ray crystallography (for complexes with human carbonic anhydrase II).     

Structural diversity in poly(disulfide)s

This article reports a synthetic methodology for single step preparation of telechelic poly(disulfide)s (PDS) by step‐growth polymerization between a di‐thiol and a commercially available monomer 2,2′‐dithiodipyridine in presence of a functional group appended pyridyl disulfide moiety as the “mono‐functional impurity” (MFI). Redox‐destructible well‐defined segmented PDSs with functional chain terminal, predicted and tunable degree of polymerization and narrow polydispersity index (<2.0) could be synthesized under a mild reaction condition. Using an appropriate MFI, PDS could be synthesized with trithiocarbonate chain terminals in a single step, which could be further used as macro chain‐transfer agent (CTA) for chain growth polymerization under RAFT mechanism producing ABA type tri‐block copolymer wherein the B block consists of the degradable PDS chain. By copolymerization between a hydrophobic di‐thiol monomer and a hydroxyl group appended di‐thiol monomer, PDS could be prepared with pendant hydroxyl functional group which was utilized to initiate ring opening polymerization of cyclic lactide monomers producing well‐defined degradable graft‐copolymer. The pendant hydroxyl groups were further utilized to anchor a polar carboxylic group to the degradable PDS backbone which under basic condition showed aqueous self‐assembly generating micelle‐like structure with hydrophobic guest encapsulation ability and glutathione responsive sustained release. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2018 , 56, 194–202     

Reverse fragment based drug discovery approach via simple estimation of fragment contributions

1. Download : Download high-res image (110KB)
2. Download : Download full-size image